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Abstract 

Photon statistics and phase properties of two-mode squeezed number states are stud- 
ied. It is shown that photon number distribution and Pegg-Barnett phase distribution for 
such states have similar (A r + l)-peak structure for nonzero value of the difference in the 
number of photons between modes. Exact analytical formulas for phase distributions based 
on different phase approaches are derived. The Pegg-Barnett phase distribution and the 
phase quasiprobability distribution associated with the Wigner function are close to each 
other, while the phase quasiprobability distribution associated with the Q function carries 
less phase information. 


1 Introduction 

Recent developments in quantum optics have led to new proposals to generate number states 
of the electromagnetic field using conditioned measurements techniques [1] or the properties of 
atom-field interactions in microwave cavities in the micromaser [2], The precisely defined two- 
mode photon number state |;V + q,N) can be used as an input field in a squeezing device, such 
as a parametric amplifier. The model involves a signal and an idler modes driven by a classical 
pump. The Hamiltonian for the two coupled modes is taken to be [3, 4] (we set k = 1) 

H = u> a a *d f uj b tib — i{gabexp(iu;t) — ^'6 t a t exp(— 

where u is the pump frequency and g is the effective intermode coupling constant . If we consider 
exact resonance u> = u a +u b then the Hamiltonian may be transformed into the interaction picture 

Hj — —2 {gab — g m & a ^) . 

In this picture the time-evolution operator is 

exp(-iHjt) = exp {—gtab + <7"^^}, 

and is immediately identifiable as a time-dependent two-mode squeezed operator: 

exp (~iHit) = S(gt). 

with squeezing parameter £ = gt . The output state at time i will be the two-mode squeezed 
number state 

10 = exp(-0/6 + r^« f )|.V + q. .V). 
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The properties of this state are phase dependent and it should be interesting to study them. 

The problem of the quantum description of the optical field phase has been the subject of 
considerable study for many years [ 5 ]. This is connected with the difficulty in constructing a 
Hermitian phase operator. Within the past few years the notion of phase variables in quantum 
systems has been greatly clarified. Pegg and Barnett [6]-[S] have shown how such an operator can 
be defined for quantized electromagnetic fields. This new formalism makes it possible to describe 
the quantum properties of optical phase in a direct way within quantum mechanics on the basis 

of the Hermitian phase operator and its eigenstates. 

A quite different approach to the concepts of the phase variable has also been widely used 
in quantum optics [9]— [1 1] and which involves quantum quasiprobability distributions such as 
the Q function and the Wigner function rather than Hermitian operators and their eigenstates. 
These quasiprobability distributions depend upon the complex eigenvalue o of the non- Hermitian 
annihilation operator, which can be expressed in terms of a radial variable \a\ and a “phase 1 9 
both of which are real. If we integrate over the radius, the resulting distributions are pei iodic 
in the phase angle and, for the most of states they satisfy all properties required by a proper 
phase distribution. In recent papers, the Pegg-Barnett phase distribution have been compared 
with those distributions obtained from the Wigner and Q functions by integrating them over the 
radius for the multi-photon down-conversion [12], displaced number states and displaced therma 
states [ 13 ], squeezed number states and squeezed thermal states [ 14 ]. In this paper we extend 
such comparison onto two-mode case. 

The purpose of this paper is to study photon statistics and phase properties of the two- 
mode squeezed number states which can be considered as a natural generalization of a two -mode 
squeezed vacuum state. 


2 Photon number statistics 

Consider two modes of the electromagnetic field, which have annihilation operators a and 
b. A two-mode squeezed number state (1MSNS) is defined by acting with the squeeze opeiatoi 
S{r,ip) on the two-mode number state | N + q,N), that is 

\N + q, AO(r.v') = S(r, ip)\N -f <7, N), q> 0, (!) 

where q is the difference in the number of photons between two modes and 

S(r,^p) = exp[r(a6e _2,i,f ’ — f^a^e 2 "^)]. (-) 

In problems in which photons are either created in pairs or destroyed in pairs the value of q remains 
constant. Note that in many applications where pair creation occurs starting from vacuum, the 
parameter q will be zero. The number state decomposition of TMSNS can be written as 

I N + q , N) { r,s) = \n + <•/, n){n + q , n\N + q , A = 

n 

= b n e'*” |n + <7, n), ( 3 ) 
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where 


and 


K = ^ h r)-^ ° V!(JV + < ' )ln!( ’ 1 + ■ ?)!) ' /I 


min(n,N) 


^ kUn - b 


(-l) n -*(sinhr)' 


- 2k 


- kl(n-k)\(N- k)\(q + k)l 


( 4 ) 


<Pn = (n- N)<p (5) 

with v? being a phase of squeezing. The above amplitude is obtained by using the factored form 
of the two-mode squeeze operator [15] 

5(r,<^) = (cosh r ) -1 exp[— a t 6 t e 2tv tanh r] 

x exp[-(a f a + b t 6) ln(cosh r)] exp[a&e~ 2,v ‘ tanh r], (6) 

The mean number of photons in the TMSNS is 

(a^a + b^ b) = (2 N + q + 1) cosh 2r — 1. ( 7 ) 

The joint probability to find n Q photons in mode a and n b photons in mode b is given bv 

P(Va,n b ) = | <77 a , 77 6 1 .V + (?,.'V)(,-.v)| 2 - (S) 

Using (3) and (4), we get 

P(n a ,n b ) = P(n + q, n ( 9 ) 

where 

P ( U + ?< n ) = F q( n ) = l^nP- (10) 

As we can see in Fig. 1, photon number distribution P q (n) has an oscillatory behaviour. Such a 
behaviour is a consequence of interference in four-dimensional phase space [16]. We would like to 
emphasize a presence of (N + 1) peaks in the photon number distribution. The similar behaviour 
of the photon number distribution was observed for the displaced number states [17]. It should 
be stressed that such a peak structure for TMSNS can be revealed only for those values of the 
parameter q greater than a certain number. This number depends on the value of N and for large 
N we ought to choose large values for such a number. Otherwise, some adjacent peaks in the 

photon number distribution might overlap and thus (N + l)-peak structure cannot be certainly 
discerned. 
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FIG. 1. Photon number distribution for the two-mode squeezed number state with r - 0.5, q 
and (a) jV = 0, (b) N = 1, (c) N = 2. 





3 Quasiprobability distributions 

In this section we examine the representation of TMSNS by quasiprobability phase-space 
distributions. For convenience we choose the squeezing parameter to be real, £ = r. The two-mode 
quasiprobability distributions are formed by a natural generalization of those for the single-mode 
fields (18]. The Glauber-Sudarshan V function, the Wigner function and the Q function are 
obtained by evaluating the Fourier transforms 

1 f 00 

= -j / d 2 t;d 2 f exp(a?/‘ - a*t>)exp (#* - (11) 

J — OO 

from the characteristic functions 

cW (?»0 = exp (f (M 2 + |£| 2 )) Tr {pexp(r ? at _ Tj m a)exp{ $ - £*6) J (12) 

where s = 1 if V = V y s = 0 if V = IF and s = -1 if V = Q. We would like to notice that 
there no exists well-defined Glauber-Sudarshan P function for states under consideration owing 
to their nonclassical nature [IS, 19]. 

According to ref. [20], the Q function can alternatively be defined as 

(13) 

From this definition we see that the Q function is always non-negative. Using the definition of 
the density matrix of TMSNS 

p = S(r)|A’ + q.N)(.\ + g.A'ISV) (14) 

and the factored form of the squeeze operator (6). we obtain 

= ^eosh J)-«* v+2 * +2 exp [“^ Q J J ‘) *anh r] exp[-(|a| 2 + |i| 2 )] 

(^inh r)"+*.V!(;V + g)!(q- j-)- v -"( Q J)- v -* 
ha n - k ''( N - n )H- V - *)!(*V + q - n)!(.V + q - k)V 

As to the Wigner function, it can also be represented as [20] 

(q^) = |p£) a (2Q)£) t (2J)exp[i-(u t a + 6 f 6)j| . (16) 

where D a (-j) and D (>(y) are the displacement operators for modes a and b respect ivelv. It is 
straightforward to evaluate the Wigner function using eq. (14) and the operator transforma- 
tions (15] 

S t (r)aS(r) = a cosh r — 6* sinh r, (17) 

i* f (r)6S(r) = tcosh r - a* sinh r. ^ 1 j;) 
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and their Hermitian conjugates. We find a quite simple analytical form for the Wigner function 

W(a,0) = — j( — 1)’ exp[— 2cosh2r(|a| 2 + \P\ 2 ) — 2sinh2r(a / 3 + a j3 )] 

xLv ((2 sinh r|o|) 2 + (2 cosh r\j3\) 2 + 2sinh2r(o5 + a’/?*)) 

xL.v+q ((2cosh r|aj) 2 -1- (2sinhr|/?|) 2 + 2sinh2r(a$ + a 3 )) . (19) 

where L n (x) is the Laguerre polynomial of order n. From eqs. (15) and (19) one can see that the 
Q function and the Wigner function depend on the sum of the phases 0 a + 6 b only. This fact clearly 
exhibits the correlated nature of the two-mode squeezed number states. In the next section we 
will employ the quasiprobability functions in consideration of phase properties of these states. 


4 Phase distributions 

Now we employ the two-mode Pegg-Barnett phase formalism [21], [22] to find the phase 
distribution function for such states. This formalism is based on the observation that the Hermitian 
phase operator can be defined in a finite-dimensional state space, spanned by the number states. 
The main idea of the Pegg-Barnett formalism is to evaluate all necessary expectation values on 
this finite-dimensional state space, and only after that the dimension of the space is allowed to 
tend to infinity. Having the number state decomposition (3) of TMSNS we can determine the 
continuous joint phase probability distribution for the continuous phase variables 0 a and 0 b . which 

is given by 

P{O a ,0 b ) = |l +2^&Acos[(n - fc)(0 a + M j < (20) 

l n>k * 


where b n are given by Eq. (4). The distribution (20) is normalized such that 



P(0 a .e b )d9 a d0 b = 1. 


(21) 


One important phase property of TMSNS is seen directly from the form of formula (20). It is 
clear that the joint probability distribution depends on the sum of the two phases only 


P(0 a J b ) = P(9+ =O a +0 b ). 


( 22 ) 


This means the strong correlations of the two modes. Integrating P(0*.0„) over one of the phases 
gives a marginal phase distribution P{0 a ) or P(0 b ) for the phases 0 3 or 0 b . which are uniformly 

distributed 

P(0J = j P(0,.0 b )d0 b = (' 23 ^ 


P(0 b ) = P(0 2 ) = 

_ it 


(24) 
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Thus the phases 9 or 9 b of the individual modes are uniformly distributed, and the only nonuni- 
formly distributed phase quantity is the phase sum 9 + = 9 a + 0 b . In Fig. 2 we plot t he Pegg- Barnett 
p ase distribution for TMSNS in polar coordinates for different values of parameter q. For nonzero 
values of q the phase distribution shows (N + l)-lobe structure, and the greater q the more dis- 
tinct lobes become. However, when q = 0 the phase distribution has only one lobe for all .V. It 
is important to notice a remarkable resemblance in a behaviour of the phase distribution and the 
photon number distribution for TMSNS: they both display the ( N + 1)- P eak structure. Another 
significant feature of the joint phase distribution is a property of the phase locking — the phase 
sum is locked to the argument of the squeezing parameter in the limit of large squeezing 121 231 
Now consider phase quasiprobability distributions which can be obtained by integrating 
quasiprobability distribution funct.ons (11) over the radial variables [9], [10], As we have noticed 
above, V function is not well-defined function for the states under consideration and therefore 
it is impossible to determine corresponding phase quasiprobability distribution. As a result of 
integration of Q{a,p) and W{a,i3) over |a| and |/?|, we arrive to the following formula: 


b h cos[{n - k)9 + ]G iV) {n,k)G {V] {n 

y ' l n>A- 




(25) 


where the coefficients G (U) (n, k) distinguish between two distributions, and they are: 
(i) for the Q function 


<?(<?)(„ M_ rK» + *)/2 + i] 

\J n\kl 


(26) 


(ii) for the Wigner function 


G (lV) (n,A :) = y^(-i) A -™2 ( l n ~* l + 2m l/ 2 


m=0 


/(m)G-m) G< °' (m ' |n - *1 + ">)• 


(27) 


rmn 7 = 1 max(n ’ *>• A " the coefficients G<*- l(n, k) are symmetrical, G< v Hn,k) = 

(t,n) and G (n,n) = 1. Note, that such expressions (20), (25) for the phase distributions 
are valid for all two-mode states with the number state decomposition like in (3). In Fig. 3 
we s ow he plots of the three phase distributions in polar coordinates for TMSNS calculated 
ccording to formulas (20) and (25) with the coefficients (26) and (27) for different values of N 
and nonzero q It is seen that the Pegg-Barnett phase distribution and P^UO.) are similar 
and have the N + 1 lobes, while P^(9 + ) is much broader and has only one lobe. In the case 
* - ° al three distributions have the same form of one lobe. So, as in the case of displaced 
p ( „ 0 b ;; * a ‘ eS I 13 !’ tl l ei i e 18 an essent,aJ difference in the phase information carried by and 

f(/+). Because of the averaging procedure with the “probabilities'’ G {Q) (n. k)G^(n + r/. k + q) 
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some phase information is lost in P^(9 + ). The Pegg-Barnett phase distribution is very dose 
to the distribution P (W) {9 + ), although it is not identical to it. The phase peaks of P { ] {0 + ) are 
slightly narrower than those of P {PB) {9+). The greater the difference in number of photons q the 
closer these two distributions. Basically they carry the same phase information. This similarity 
is in agreement with the area of overlap in phase space arguments, which are that the Wigner 
function represents quantum states in the phase space [10]. However, the Wigner function can 
take on negative values and the positive definiteness of P( w \0 + ) is not automatically guaranteed, 
while there are no such problems with the Pegg-Barnett phase distribution. 


5 Conclusions 

We have discussed photon statistics and phase properties of the two-mode squeezed num- 
ber states showing that the photon number distribution and the Pegg-Barnett phase distiibution 
for such states exhibit the similar N + 1-peak structure for nonzero values of the diffeience in 
the number of photons q between modes. We have compared the Pegg-Barnett phase distribu- 
tion with the phase quasiprobability distributions P (Q ' > (0+) and P (lv, (#+) obtained by integrating 
the Q function and the Wigner function over the radial coordinates. We have shown that the 
Pegg-Barnett phase distribution and the distribution P^ w \9+) carry basically the same phase 
information, while the distribution P (Q) {9 + ) loses an essential part of the phase information. 
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